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Abstract 

In this paper we examine the zero and first order eigenvalue fluctuations for the 
/3-Hermite and /3-Laguerre ensembles, using the matrix models we described in [S], 
in the limit as /3 — > oo. We find that the fluctuations are described by Gaussians 
of variance 0(l//3), centered at the roots of a corresponding Hermite (Laguerre) 
polynomial. We also show that the approximation is very good, even for small values 
of /3, by plotting exact level densities versus sum of Gaussians approximations. 



1 Introduction 



This paper provides insight into the shape of random matrix laws such as the finite 
semi-circle law, the finite quarter-circle law and its generalization. We begin with 
a simple example. Suppose A is a random k x k complex matrix with real and 
imaginary parts all i.i.d. standard normals. Let S = {A + A^)/2 be the Hermitian 
part of A. The matrix S has a distribution commonly known as the Gaussian 
Unitary Ensemble; this matrix distribution and the joint distribution of its (real) 
eigenvalues have been well studied. For a good reference on the subject, see Mehta 

nni- 

We draw below histograms of normalized eigenvalues taken from this distribu- 
tion, the known theoretical distribution (see |1(J[ page 93]), and the semicircle limit 
corresponding to fc — > cx). For the histograms, we have chosen 40, 000 samples from 
the GUE with fc = 4 and fc = 6. 



k=4, p = 2 



k=6, 




Figure 1: Histograms of eigenvalues, finite {k = 4, 6) exact level densities, and semicircle 
law (k — > oo) for matrices of the GUE ensemble 



Notice, for k finite, the k "bumps" in the distribution that wiggle above and 
below the semi-circle. A natural question to many engineers, physicists, mathe- 
maticians, and other scientists who have seen these pictures is whether they can 
be well approximated by the sum of k appropriately chosen Gaussians. (Of course 
when fc = 1, this is exactly true.) The answer, as proved in this paper, is yes. We 
give a sum of Gaussians approximation that is asymptotically correct for the /3 — )■ oo 
limit but useful even for small values of (3. 

For those well versed in random matrix theory, the GUE is the (3 — 2 case 
of a Hermite matrix ensemble 10 . Had we started with A real (quaternion), we 
would have the Gaussian Orthogonal (Symplectic) Ensemble corresponding to /? = 1 
(^ = 4). 

The joint eigenvalue density f^{Xi, . . . ,Xk) defined on R*-' for the k eigenvalues 
for an arbitrary /3 > is given in the formula below. 

/,^(Ai,...,A.) = (2.)-'=/^n-ni±|^ n |A.-A,|^e"-'^''. (1) 

j = l ^ U + .? 2 l<i<j<k 

Note the "repulsion" factor A(A) ee A(Ai, . . . , Afc) = Jl |Aj - Aj|. 

l<i<j<k 

Similarly, for the kxk Laguerre ensembles of statistics (Wishart matrix theory), 
the joint eigenvalues density ^ is defined on [0, oo)'^ for arbitrary /3 and parameter 
a > {k — l)/3/2, and is given below; once again note the repulsion factor A(A): 

fi,dXu...,Xk)=ci^ n \K->^A'tlK-''-'''^'-'e-S^'''\ (2) 

l<i<j<k i=l 



where 



C/3,, 



-ka 



r(i + f) 



In [S] , we have found (real) tridiagonal matrix models whose eigenvalue distribu- 
tions are given by and we depict the distributions in Tabled Note that the 
variables have either standard normal distribution or a x distribution (sometimes 
scaled by -\/2). 

For generating efficiently eigenvalues for the /3-ensemble distributions, we rec- 
ommend using the tridiagonal/bidiagonal model above. 

The marginal density of a single eigenvalue (known as the level density) can be 
computed in the case of the Hermite ensembles for f3 an even integer |2] using the 
MOPS software [H]. For general /3 perhaps a little more research is needed, but it 
is likely that computational techniques are not far out of reach. 

For fixed k and general /? one finds that for /3 getting larger, the bumps of 
Figure ^ get "bumpier" . To be precise, for the Hermite ensembles, we prove here 
that at /3 = oo the bumps become delta functions at the roots of the kth Hermite 
polynomial, while for (3 large, the bumps behave like Gaussians centered at these 
roots with variance 0(l//3). 

The model of (3 as an inverse temperature is apparent from As (3 goes 

to 0, the strength of the repulsion factor A(A) decreases until annihilation; the 
interdependence among eigenvalues disappears, and the randomness increases (each 
eigenvalue behaves like an independent Gaussian). In the frozen state (/? = oo), 
we can imagine the k eigenvalues fixed at the roots of the Hermite polynomial. 
Warming the system a little (/3 very large but not infinite) gives the particles a 
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Hermite matrix 
n G N 


1 


' iV(0,2) X{n-X)p \ 
X(n-l)/3 ^(0,2) X(n-2)/3 

X2/3 iV(0,2) X/3 
, X/3 iV(0,2) / 


Laguerre matrix 

m G N 
a G M 
a > |(m - 1) 


Bp- 


g , where 

X/3(m-l) X2a-/3 
\ Xfi X2a-/3(m-l) / 



Table 1: Tridiagonal matrix models for the /3-Hermite and /3-Laguerre ensembles with 
any /3 > 0. 



little energy, and the eigenvalues have Gaussian distribution to first order around 
the Hermite polynomial roots. 

Similarly, in the Laguerre case, at /3 = 0, the eigenvalues become i.i.d. variables 
with distribution xia- P grows the eigenvalues have Gaussian distribution to 
first order around the Laguerre polynomial roots, while aX (3 = oo we reach the 
freezing point when the eigenvalues are fixed at those roots. 

In the following, we use first order eigenvalue perturbation theory and the tridi- 
agonal ensembles in |S] to rigorously investigate this phenomenon mathematically 
obtaining precisely the asymptotic variance along with the mean. 

These results draw a parallel to the Tracy- Widom laws |13[I14| for the [3 = 1, 2, 4- 
Hermite ensembles, later extended to /3 = 1, 2-Laguerre ensembles by Johansson |S] 
and Johnstone 

The Tracy- Widom laws compute the fluctuation in the distribution of the largest 
eigenvalue of a /3-Hermite ensemble with /3 = 1, 2, 4, as fc ^ oo, and obtain it in 
terms of the solution to a Painleve differential equation. From the semicircle law, 
we know that as n ^ oo, regardless of (3, the largest eigenvalue (scaled by \/2kP) 
goes to 1. From Theorem 13.11 constrained to i — 1, with the help of (TTI and ^ 
page 450], we obtain Corollary 13.31 which gives an intuition of how the (3 = 1,2,4 
Tracy- Widom distributions evolve towards a normal distribution at (3 — oo. 

The theoretical results of Section|31are similar to the "Central Limit Theorems" , 
i.e. the computation of the global fluctuations from the semicircle and semicircle- 
type laws done by Johansson in 7 for Hermite-like ensembles of any (3 and by 
Silverstein and Bai 12 for a class of Laguerre-like ensembles with real or complex 
entries {(3 = 1,2). Roughly said, the eigenvalues can be thought of as fluctuating 
(like Gaussians) around the roots of the corresponding orthogonal polynomial as (3 
grows large; if one lets n grow large, too, the global eigenvalue fluctuation becomes 
a Gaussian process. The larger 1//3, the "warmer" it gets, and the larger the 
"vibration". The larger /3, the "cooler" it gets, and the eigenvalues "freeze" into 
place. 



3 



At the end of Section Owe perform computational experiments to see how good 
the j3 large approximation is even for relatively small k and (3. 



2 Eigenvalue perturbation and x asymptotics 

In this section we present two lemmas we need in the proofs of our main results 
(Theorems O and ETJ . 

The first lemma involves perturbation theory; for a good reference on Pertur- 
bation theory and a more general form of the result below, see Demmel's book 
Section 4.3]. 

Lemma 2.1. Let A and B be n x n symmetric matrices, and let e > 0. Assume A 
has all distinct eigenvalues. Let M = A-\-eB-\-o{e), where by o(e) we mean a matrix 
in which every entry goes to faster than e. Let Xi{X) denote the ith eigenvalue of 
X , for 1 < i < n. Finally, let Q be an eigenvector matrix for A. Then 

lim i(A,(M) - A,(A)) = Q(; , ifBQ{:, i) , 

where, following MATLAB notation, Q{:,i) represents the ith column of Q. 
Remark 2.2. Equivalently, for every 1 < i < n, 

X,{M) = K{A) + eQ(; , ifBQ{:, i) + o(e) . 

The second result is an approximation lemma for the Xr distribution as r grows 
large. 

Lemma 2.3. Let r > 0, and let X be a variable with distribution Xr- Then as 
r ^ oo the distribution of X ~ ^Jr converges to a normal distribution of mean and 
variance 1/2. 

Proof. We prove this lemma by looking at the density function of Xr when r ^ oo. 
First we will show that E[X] — %/r ^ as r ^ oo. Since 

2i-iV2 v{^) r(^) 
r(§) 2i-('-+i)/2 = ^T{f^ • 

Using the Stirling approximation formula 

r(z) ^ ^-i/2e-V2^(^l + -L + o(^)^ , (4) 

for T large, Q becomes 

E{X\ = (l + 0(r-i)) . 

So as r ^ 00, the mean of X — ^fr approaches 0. Now let us examine the p.d.f. 
oiX-^. Denote by F = X - ; the p.d.f. of Y is 

We examine this p.d.f. in a "small" neighborhood of 0, such that t = o(r^/^). 
Again, we use the Stirling approximation for the Gamma function at infinity, to 
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obtain 



2l-r/2 ^r/2-1 ^ ly-'^ ^-t^ /2-^t-r/2 



1 A. t 



/TT 



and by using the Taylor series for (r — 1) ln(l + t/^/r) and the fact that t = o(r^/^), 
we obtain 

m = ^e-^' (l + o(±]] . (5) 



Thus, on any fixed interval, the p.d.f. of X — ^/r converges to the p.d.f. of 
a centered normal of variance 1/2. This is enough to prove that the c.d.f. of 
Y = X — y/r converges to the c.d.f. of a centered normal of variance 1/2. □ 



3 /3-Hermite: zero and first-order approximations 



, hj^' be the roots of the fcth imivariate Hermite 



Let k be fixed, and let ft,^'^-' , . . . , h^j^^ 
polynomial (where the Hermite polynomials Hq{x), Hi(x), . . . are orthonormal 
with respect to the weight on (—00,00)). 

Let A/j be a random matrix from the /3-Hermite ensemble of size k, scaled by 
l/v^Sfc^. For the remainder of this section, we think of /? as a parameter. 

We state and prove the following theorem. 

Theorem 3.1. Let Xi{A/}) be the ith largest eigenvalue 0} Ap, for any fixed 1 <i < 
k. Then, as /3 —> 00, 



KiAff) 



2k 



(fe) 



Moreover, as /3 — > 00, 



^^^\...,Xk{A0)- 



2k 



2k 



G 



where G = (Gi, G2, • • • , Gfc) is a k-variate Gaussian with covariance matrix 



Cov{Gi, Gj) = 



T Hl{hf^)Hl{hf) + '^ Hi^,{hf^)Hi{hf^)m+,{^^^^ 
(=0 ;=o 



/fc-i 



\l=0 J \l=0 



k-1 



Proof. Let H be the k x k symmetric tridiagonal matrix 



H 



( s/k^ 

y/k^ y/k^ 

\fk^ 



72 



\/T 
sfl / 



(6) 
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This matrix is the tridiagonal matrix corresponding to the 3-term recurrence 
for Hermite polynomials. Its eigenvalues are the roots of the fcth Hermite poly- 



nomial Hk [x) (recall that we denoted them by /i^'^-' , . . . , h^^' ) , and the eigenvector 
corresponding to the ith eigenvalue h^^^ is 



Hk-2{h 



V Hoihl"^) ) 

Lemma 3.2. Let Ap be as defined in the beginning of this section. Then almost 

surely 

lim V2fc/3A/3 - = Z , 

/3— >oo 

where Z is a tridiagonal matrix with standard normal variables on the diagonal 
and normal variables of mean and variance 1/4 on the suhdiagonal. All normal 
variables in Z are mutually independent, subject only to the symmetry. 

From now on we use the notation 



Z 



Nk-i 



V 



Nk-1 
Mk-i 

Nk-2 



Nk-2 
Mk-2 



M2 Ni 
Ni Ml ) 



(7) 



with Mi being independent standard normals, while Ni are independent normals 
of mean and variance 1/4; all normal variables in Z are mutually independent, 
subject only to the symmetry. 

Lemma 13.21 follows immediately from Lemma 12.31 since we are dealing with a 
finite number (fc — 1) of x variables on the sub-diagonal of A^, each converging in 
distributions to a normal variable. 

Hence we have that, entry by entry, 



Af- 



^=Z + o( , , 



in distributions, as /3 — > 00. 

Thus all zero- and first-order properties of Ap are the same as for the random 
matrix -^^H + --^==pZ, where Z as as above. In particular, for any I < i < k, 



Xi{Ap) = A, 



H+^^Z] +o( ^ , , 
2k V2kp ) \VPJ 



in distributions as /3 — > 00. 

Finally, with the help of perturbation theory Lemma |2. II we obtain that for any 
l<i<k, 



2k 



/2kp vfv. 



in distributions as /3 — > 00. 
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Hence, using the notation ((Tj), 

The statement of Theorem 13.11 fohows . □ 

Letting fc — > cx) in Theorem 13. II we obtain the CoroUary below. 

Corollary 3.3. Let Ap be a matrix from the k x k f3-Hermite ensemble, scaled by 
^2k/3 ' '^^'^ -^1(^/3) ^6 tf^^ largest eigenvalue of Ap. Then 

lim Um fc"2/3 (Ai(A/3) - 1) ^ ^ + a^G, 
where Oi = —2.38810 . . . zs t/ie largest root of the Airy Ai function (see 'Jj), and 

i4 



a^^2^ ^ ^ - 0.41050. 



ki?{x + ao) dx) 



2 



Proof. The corohary follows by using the special properties of the Hermite polyno- 
mial roots and the Airy function, as in ^ and |11[ (1.81), (6.32)], and from the fact 
that the matrix lO is a discretization of the Airy operator, and as fc ^ cx), the first 
eigenvector of the matrix becomes the Ai solution of the operator which is bounded 
at cxD and positive. □ 

4 /?-Laguerre: zero and first-order approximations 

ik') (k) 

Let k be fixed. Given a fixed 7 > 0, let ll ,...,lj, be the roots of the fcth 
Laguerre polynomial of parameter 7 — 1, ^. Note that for any 7 > 0, the 
Laguerre polynomials Lq, LJ, . . . are orthonormal with respect to the weight x'^e^^ 
on [0,00); for 7 S [—1,0) they admit formal definitions. 

Let 3/3 be a random matrix from the /3-Laguerre ensemble of size k and pa- 
rameter a/3, scaled by l/fc/3. For the remainder of this section, we think of /3 as a 
parameter. Suppose that, as /3 grows large, 

hm ^^i(fc + 7-l). 

Note that the requirement ap>{k— l)/3/2 constrains 7 to be positive. 

Theorem 4.1. Let \i{Bfj) be the ith largest eigenvalue of Bp, for any fixed 1 < 
i < k. Then, as (3 00, 



Moreover, as P —f 00, 



^ (^Ai(B^) - ^l[''\ X2{Bp) - ^lf\. . . , Xk{Bp) - ^4 

where G = (Gi, G2, . . . , Gk) is a centered k-variate Gaussian of covariance matrix 

h + k- l)[Ll_Sf^)f{Ll_,[lf^)f + A^{i,j) + B,{z,j) + Gk{i,j) + D,{i,j) 



Cov{G^, Gj) = 2- 



1=0 J \i=a 



7 



fe-1 



5^(7 + 2(fc - - l){Ll_,_Sf^)fLU_Sf )f , 
1=1 



1=1 



and 



Proof. The proof follows in the footsteps of that of Theorem 13. II 
Let L-y be the k x k (symmetric) positive definite matrix 



/ 7+fc- 1 



V7 + fc - IVfc^ 



V7 + fc- l^fc^ 2(fc - 2) + 7 + 1 + - 2y/k^ 

V-y + k- 2y/k^ 2(fc-3)+7+l 



V7T2./2 

^/7T2^/2 3 + 7 V7n\/T 

vtttv/t 1+7 y 



(8) 



We can write = B^B^ , with 

/ V7 + fc ~ 1 



V 



72 V^TT 



(9) 



The matrix has as eigenvalues the roots of the A;th Laguerre poly- 
nomial of parameter 7 — 1, L^~^(x) (recall that we have denoted them by 
. . . , s-'^d an eigenvector corresponding to the ith eigenvalue I'f'^ is 

Wi = : 

L?(/f)) 

We define 0i = tfj/l I2 to be a length 1 eigenvector corresponding to 
the ith eigenvalue U. 
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Lemma 4.2. Let Bf^ he as in the statement of Theorem \4.1\ Then 



1 

V2' 



lim kf3Bf^ - PL^ = ^{B^Z^ + ZBI^) 



in distributions, where Z is a lower bidiagonal matrix with standard normal 
variables on the diagonal and on the suhdiagonal. All normal variables in Z 
are mutually independent, subject only to the symmetry constraint. 

We use the notation 

\ 



Nk^i Mk-i 



V 



N2 M2 

Ni A h I 



(10) 



Once again, the proof for this lemma follows from the construction of the 
Laguerre matrix as a lower bidiagonal random matrix times its transpose, 
and from Lemma 12.31 applied to the x entries on the bidiagonal random 
matrix (there is a finite number 2/c — 1 of them). 

Just as in the Hermite case. Lemma l4 . 2 1 allows us to write that, entry by 
entry, 



L^ + 



1 



{B^Z^ + ZB') + o 



1 



{'VP J ' 



in distributions, as /? — > cxo. Thus once again. 



1 



1 



X,{Bp) = Xi{-L^ + j^iB.Z-^ + ZB'^)) + 



in distributions, as /3 ^ 00. 

Finally, perturbation theory Lemma l2 . II applies once again to yield that, 
as /3 — > 00, 



lAk)^ 1 wJ{B,Z^ + ZBi;))w, 



wjwi 



+ 



in distributions. 

,JB, 



Since wfB-.Z'^Wi = wf ZB'^Wi, we can write that, as /3 ^ 00, 



1 



V2 wjB^Z^ 



^'^^'^ - k'^ ' kV]3 wfw, 

in distributions. 

Thus, using notation we can write 



Wi 



+ 



1 



[k) 



V2 ^{Ll{lf^)f + Sumi + Sums 



z=o 
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with 



fc-i 



Sumi = j;(y7T7(L7(/f))2 + V/L/(/f))L7_,(/f)))M,+i, and 
/=i 
fc— 1 

Sum2 = E(V7T7L7(/f)L7_,(/f) + V/(L7_,(/f)))^)iV,. 
1=1 

The statement of the theorem fohows. 



□ 



5 Applications: Level densities 

We can compare the large P asymptotics to the theoretical answer for the distribu- 
tion of a randomly chosen eigenvalue. For large n, this is the well-know semicircle 
law (for the Hermite ensembles) or equivalent thereof (for Laguerre ensembles) , but 
we are interested in finite n. 

We found that even for /? small, the approximation can be quite reasonable. 

We summarize the large /3 answer as a sum of Gaussians in Corollaries 15 . II and 

o 

Corollary 5.1. Letk be fixed, and fk,p be the level density of the scaled (by l/^/2kj3) 
k X k (3-Hermite ensemble. Let gk^p be as below: 

1 . ^ ^ 1 (J^-Mj)^ 



k ^ V2^cr, 



^(fc) 

where = and Ui — ,^ \/ Var{Gi), with hi and Var{Gi) as in Section\^ 

Then for any x, 

lim \/P {fk.p{x) - gk.p{x)) = 0. 

Corollary 5.2. Let k and 7 > fee fixed, and fk.fj,-y be the level density of the scaled 
(by l/(fc/3)j k X k f3-Laguerre ensemble of parameter a — ^(fc — 1 + 7). Let gk,i3,'y 
be as below: 

''^'-^'^ ^ k^^V^.' 

where iii = and at = -^j^\J yo^riGi), with l^^'^ and Var{Gi) as in Section^ 
Then for any x, 

lim ^ {fk,0^^{x) - gk,f3.-({x)) = . 

While these approximations are simple enough (a sum of Gaussians is an easily 
recognizable shape that is also easy to work with), one may wonder how big /? has 
to be in order for these approximations to become "accurate" (for example, in order 
to appear accurate in a plot, the approximations have to be accurate to about 2-3 
digits). We have found that, in either of the two cases, the answer is surprisingly 
low. 

In the following two subsections, we have used only even integer values of /3 for 
our plots, because (in addition to (3 — 1) those are the only ones for which (to the 
best of our knowledge) there are exact formulas for the level densities. The plots 
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were obtained with the help of our Maple Library, MOPs {Multivariate Orthogonal 
Polynomials (symbolically)), which was used for computing the orthogonal and Jack 
polynomial quantities involved; these were translated into polynomials which were 
then plotted in MATLAB. For a reference on MOPs see 

5.1 Level density plots: the Hermite case 

In the following, we illustrate the accuracy of the level density approximation by a 
sum of Gaussians for f3 relatively small (4 to 10) by plotting them against the true 
level densities. 

Figure El plots the Hermite case with fc = 4. 



k = 4, |3 = 6 




-1.5 -1 -0.5 0.5 1 1.5 -1.5 -1 -0.5 0.5 1 1.5 



k = 4, (3= 10 




-1.5 -1 -0.5 0.5 1 1.5 -1.5 -1 -0.5 0.5 1 1.5 



Figure 2: Hermite case: sum of Gaussians approximation to the level densities (dots) 
and exact level densities (lines) for = 4, and /3 = 4, 6, 8, 10 

In Figure 121 we let k — 4, and gradually increase (3 (from 4 to 10) to show how 
the plots become more and more similar. For /3 = 10, the two plots appear to 
overlap. 

If we plot the densities for k = 7 (as in Figure O, (3 = 6 already provides a very 
good approximation. 

We can conclude that the approximation works well for low values of /3, in the 
Hermite case. 

5.2 Level densities: the Laguerre case 

In the Laguerre case, we cut the parameter cube with two different slices, as ex- 
plained below. For plotting purposes we have considered /c = 4 in both. 
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k = 7,p = 4 




k = 7, |3 = 6 




Figure 3: Hermite case: sum of Gaussians approximation to the level densities (dots) 
and exact level densities (lines) for k = 7, and /3 = 2, 4, 6. 



To specify a Laguerre density, there is always an intrinsic parameter: the power 
"p" (as in xPe~'^). However, in this story, there are two Laguerre densities: one 
for the eigenvalue p.d.f. that is in the Laguerre ensemble density, and a second 
(different!) one for the Laguerre polynomial corresponding to the limiting level 
density. We found this surprising at first. In order to study /3 — s- cx) limits we 
looked at both possibilities of holding a parameter constant, as depicted in the 
table below. 





Fixed 
quantities 


Variable 


Other 
quantities 


Eigenvalue 
p.d.f 


Limiting 
Laguerre 
polynomial 


a) 


k, 7 


/3 ^ oo 


a= |(fc + 7-l) 






b) 


k, P 


/3 -> cx) 


a=p+f(fc-l) 
7= §(P+1) 


ciAi^nliAre-^-/2 





Case a). This case holds 7 (and therefore the limiting Laguerre polynomial, whose 
roots are the limits of the scaled eigenvalues) constant as /3 ^ 00. 

Note that both the Laguerre ensemble parameter a = '|(fc + 7— 1) and the power 
j3 = 7^ — 1 are increasing functions of (3. 

By prescribing 7, in the limit as /3 — > 00 the plot should become a sum of delta 
functions at the roots of the Laguerre polynomial L^^^(a;). 
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In Figure H we take fc = 4, 7 = 1, /3 = 4, 6, 8, 10, and a = 8, 12, 16, 20 (equiva- 
lently, p = 1, 2, 3, 4). Note that the approximation is very good for /3 = 10. 



k = 4,[3 = 6,y= 1, p = 2, a = 12 
k = 4,p = 4, Y=1,p = 1,a = 8 ^ 




0.5 1 1.5 2 2.5 3 3.5 4 0.5 1 1.5 2 2.5 3 3.5 4 



Figure 4: Laguerre case a): sum of Gaussians approximation to the level densities (dots) 
and exact level densities (lines) for A; = 4, 7 = 1, and /? = 4, 6, 8, 10 



Case b). This case holds the power p constant in the weight |A(A)|'^ IliLi ^^/'^^ 
thereby changing the parameter 7 and the Laguerre polynomial. In this second test, 
as /3 ^ 00, 7 = |(p+ 1) ^0. 

Thus as ^ 00, the plot should become a sum of delta functions at the roots 
of the polynomial L~^{x). 

The approximation works, once again, surprisingly well, as demonstrated by 
Figure where n = 4, p = 1, /3 = 4,6,8,10, and 7 = 1,2/3,1/2,2/5 (or a = 
8,11,14,17). 

Remark 5.3. Note that in this case, the smallest root converges to (which is 
the smallest root of the Laguerre polynomial L^^{x) ), and the presence of the delta 
function at in the sum of Gaussians is very clearly visible. 

Thus we can conclude that in both cases, a good approximation is obtained even 
for /3 relatively small. 
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k = 4, (3 = 6, p = 1,Y=2/3, a = 11 



k = 4, p=4,p = 1,Y =1,a = 8 




0.5 1 1.5 2 2.5 3 3.5 4 



k = 4,p = 8,p = 1,Y=1/2,a = 14 




0.5 1 1.5 2 2.5 3 3.5 4 

k = 4, p = 10, p= 1, y= 2/5, a= 17 




0.5 1 1.5 2 2.5 3 3.5 4 



0.5 1 1.5 2 2.5 3 3.5 4 



Figure 5: Laguerre case b): sum of Gaussians approximation to the level densities (dots) 
and exact level densities (lines) for = 4, p = 1, and /3 = 4, 6, 8, 10 
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6 Circular Ensembles 



Similar to the /3-Hermite ensemble, we have the circular ensembles defined by the 
joint eigenvalue e**^ (with Oj £ [0, 1]") density proportional to 

l<j<l<k 

The (3 = 2 circular ensemble is also known as the Haar measure on the unitary 
group Un- The eigenvalues of [/„ appear to be almost uniformly distributed on the 
unit circle (see the experiment with fc = 100 in Diaconis' paper 0]). For any fixed 
fc, as /3 ^ CX3, the eigenvalues freeze into place uniformly at the fcth roots of unity. 
We believe that the same Gaussian phenomenon will hold, and the fluctuation of 
eigenvalue i will behave like a normal centered at the zth root of the unity, with 
variance depending on ^. 
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